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We show that Vassiliev invariants of knots, appropriately generalized to the spin network context, 
are loop differentiable in spite of being diffeomorphism invariant. This opens the possibility of 
defining rigorously the constraints of quantum gravity as geometrical operators acting on the space 
of Vassiliev invariants of spin nets. We show how to explicitly realize the diffeomorphism constraint 
qq , on this space and present proposals for the construction of Hamiltonian constraints. 
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I. INTRODUCTION 



Quantum gravity presents a fundamental problem: it is a theory with an infinite number of degrees of freedom, 
and yet it is "topological" in the sense of being invariant under diffcomorphisms. This makes the description of the 
theory particularly troublesome. In particular, if one recourses to structures of intrinsic topological nature, they 
tend to have a discrete character that makes them incompatible with the demands of a field theory with an infinite 
■ number of degrees of freedom. To put this discussion in a more concrete setting, let us concentrate on the canonical 
formulation of quantum gravity in terms of the new variables introduced by Ashtekar. In terms of these variables 
one can construct a representation of the quantum theory purely in terms of loops (the loop representation) , and a 
suitable basis for describing quantum states is given by spin networks. That is, quantum states are labelled by spin 
networks, multivalent graphs embedded in three dimensional space, with a system of weights associated to each side 
of the graph. If one imposes the diffeomorphism constraint, one considers functions of the diffeomorphism class of a 
spin network. There exists a precise sense in which one can endow such a space with an inner product 0, in terms 
of which spin network states are orthonormal. 

In this context, there exists a proposal for the action of the Hamiltonian constraint of quantum gravity due to 
Thiemann ||. In such a proposal the Hamiltonian constraint weighed by a lapse is just given by a topological operator 
acting at the vertices of the spin networks. The Hamiltonian constraint commutes with itself, as it is expected it should 
happen in a diffeomorphism invariant context. The construction also allows for the inclusion of matter in the theory in 
a well defined regularized framework. It has been observed ||, however, that if one considers Thiemann's Hamiltonian 
as acting on a larger space of functions (that are not diffeomorphism invariant), the Hamiltonian constraint continues 
to commute with itself. Although one can arrange the right hand side of the commutator to also vanish, it appears 
that the price to pay is tantamount to having a degenerate metric B . It also appears that this result is rather generic, 
holding for many possible detailed forms of the action of the Hamiltonian at vertices. Other non-local proposals also 
seem to suffer of difficulties with the constraint algebra || . 

On the other hand over the last few years, a variety of formal results have been obtained on a different space of 
states, that in which the loop derivative is well defined. In particular, a proposal for the Hamiltonian constraint of 
quantum gravity in terms of loop derivatives exists, and it has been shown that at a formal level the classical Poisson 
brackets are reproduced by the quantum theory J(|. The main drawback of this space of states is twofold: on the one 
hand there is the fact that the loop derivative does not appear to exist on generic diffeomorphism invariant states. 
This is due to the fact that such states change discontinuously when one changes the loops and therefore one cannot 
introduce a differential operator in loop space. Moreover, definitions of the Hamiltonian in terms of the loop derivative 
have only been attempted in the context of multiloops (not spin networks). In this context, each type of loop and 
intersection has to be treated individually, and therefore most results (for instance proofs that certain states were 
annihilated by the constraints) were only of a partial nature. 

The purpose of this paper is to notice that there exist a set of especially important knot invariants that appear to 
have the property of being loop differentiable. As a consequence, one can operate on them in terms of the Hamiltonian 
and diffeomorphism constraints of quantum gravity written in terms of loop derivatives in a more systematic way. 
The invariants in question are the Vassiliev invariants, which in turn are conjectured Q to be complete enough to 



'Associate member of ICTP. 



1 



be able to separate knots. Moreover, they have recently been generalized to the spin network context ||. We will 
show explicitly that these invariants are annihilated by the diffeomorphism constraint of quantum gravity written in 
terms of loop derivatives. We will also start the analysis of the action of the regularized Hamiltonian constraint of 
the theory. All of this will be done in terms of spin networks, which allows us to discuss in a unified way all types of 
intersections and loops. 

The organization of this paper is as follows. In section II we discuss the Vassiliev invariants and their loop 
differentiability. In section III we discuss the action of the diffeomorphism constraint. In section IV we analyze 
proposals for the Hamiltonian constraint. 



II. PRELIMINARIES 



A. Quantum states 



The kind of states we will consider as candidates of interest for physical states in quantum gravity has a long 
motivation in attempts to try to show that there is a loop space representation counterpart of the Chern-Simons 
state in quantum gravity (see (9) for references). To briefly summarize, one considers general relativity written in 
terms of Ashtckar's new variables, which consist of a canonical pair given by a set of densitized triads Ef and a 
(complex) SU(2) connection A l a . If one considers a quantum representation in which wavefunctions are functions of 
the connection it is known that if one considers a wavefunction given by, 

®[A] = exp f^- J Tr(A A dA + ^A AAAij = exp (ikScs) (1) 

is annihilated jn^,[ll) by all the constraints of quantum gravity with a cosmological constant A = 24-rr/k. 

One can now consider a loop representation, obtained by expanding the states of the connection representation in 
terms of a basis of states given the traces of the holonomies built with the connection (Wilson loops). In terms of 
such a basis, the wavefunction we considered above would be given by, 

*[ 7 ] = f DAexp(ikS C s)W(rf,A) (2) 

where 7 is a loop. This integral can be thought of in a different context, as the expectation value of the Wilson 
loop in a Chern-Simons theory. Such integrals have been the subject of a lot of studies. It is known that the 
integral is related with a knot polynomial that is a regular isotopy invariant (a framing dependent knot invariant) 
known as the Kauffman bracket. It is also known that all the framing information is concentrated in an overall factor 
proportional to the exponential of the self- linking number of the knot. If one removes this overall factor one is left with 
a polynomial that is closely related to the Jones polynomial, which is ambient isotopic invariant (a genuine invariant 
under diffeomorphisms) pjfl . These polynomials appear evaluated for a particular value of their variable, usually 
denoted as q. The particular value is q = exp(2ni / k) , with k being the coupling constant of the Chern-Simons theory. 
The resulting expressions are presented not as polynomials but as infinite series in 1/fc. The coefficients of these 
series are known to be Vassiliev invariants [Q,[l3 14 1. Although these results have never been derived using rigorous 
measure theory, they have been used extensively to construct solutions to the constraint equations. It was noted that 
one could write the constraints of quantum gravity in the loop representation in terms of differential operators in 
loop space involving the loop derivative, and that one can (formally) act with these operators on states as the one 
listed above. The end result of these manipulations was that one can show that the exponential of the self-linking 
number is annihilated by the Hamiltonian constraint with a cosmological constant and that the second coefficient of 
the expansion in of the Jones polynomial is annihilated by the vacuum (zero cosmological constant) Hamiltonian 
constraint. This coefficient is a Vassiliev invariant. 

All these results had several shortcomings. When one works in the loop representation one needs to ensure that 
wavefunctions satisfy a series of identities stemming from the fact that the loop basis consists of traces of SU(2) 
matrices. These identities are called Mandelstam identities and in general imply very complex relations between 
the values of functions of loops. Therefore if one wants to consider a knot invariant as a possible candidate for a 
state of quantum gravity in the loop representation, one needs to ensure that such invariant satisfy the Mandelstam 
identities. Although this was the case for the second coefficient we discussed above, it appeared as unlikely that 
higher coefficients would satisfy the identities. Moreover, the proofs that a certain coefficient was annihilated by the 
constraint were performed for very special types of intersections (the Hamiltonian constraint only acts nontrivially 
at intersections) . Part of the intention of this paper is to show that similar results are present in the spin network 
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context, where all types of intersections can be treated in a unified fashion. A key to generalizing these results has 
been the recent observation that ambient isotopic invariants (genuine invariants under diffeomorphisms) can also be 
defined in terms of spin networks pi. 



B. Spin networks 



Spin networks are constructed considering graphs that are embedded in three dimensions. The graphs are multi- 
connected with intersections that can be tri-valent or of higher valence. Each connecting line is associated with a 
holonomy of the connection A l a in a given representation of the group in question (in our case, SU(2), representations 
are labelled by a (half)integer). One can construct a generalization of the trace of the holonomy (Wilson loop) 
which we call Wilson net, by considering the traces of the holonomies along the Wilson nets joined by appropriate 
"intertwiners" at the intersections. The intertwiners consist of invariant tensors in the group. There are several 
possibilities of how to connect the invariant tensors and the holonomies and this has led to different conventions in the 
definition of the spin networks. The convention we will choose will follow closely that of Witten and Martin [ ]l5| , |l6| |, 
and differs of those of other authors in the field Jl7],||,[l8| . 

Consider a multi-connected graph embedded in three dimensions. In order to associate a holonomy to each line in 
the graph, we need to associate an orientation to each line, so from now on we will consider oriented graphs, with each 
line associated with a number which labels the representation in which we are considering the holonomy along such 
a line. The next task is to assign the intertwiners at each intersection. Let us consider trivalent graphs. For graphs 
of higher valence, one can choose a decomposition of the higher order intersections into superpositions of graphs with 
trivalent ones p7[ . At each intersection one can have several possibilities, depending on if the lines are all incoming, 
all outgoing or some are ingoing and some outgoing. Let us consider the case of either all ingoing or all outgoing lines. 
In those cases, there exist two ways of contracting the holonomies with the 3jm symbol which intertwines them, given 
that the latter is cyclic in its three entries. We will label one of the possibilities "+" and the other "— " and we will 
refer to them as the "orientations" of the vertex. For a positive orientation we therefore have that the intertwiner is, 




Jl 32 33 
mi 7712 7773 



(3) 



whereas for a negative orientation we have, 




3i 33 32 

mi 7713 7772 



(-1) 



H+h+h I -?2 33 

mi 7712 77T-3 



(4) 



It is useful to introduce a tensor notation. In this notation we denote the intertwiners as V and use lower indices to 
denote outgoing lines whereas upper indices denote ingoing lines. Therefore, 



T/ _ T/Jimi J2m 2 j 3 m 3 _ 

V Jimi J2TO2 ' 



3i 32 33 
mi m-i 7773 



(5) 



Holonomies are also associated with tensorial objects, which we denote by U. These objects are labelled by the 
appropriate representation of the group j and have two indices, the upper index is associated to the point of departure 
and the lower index to the point of arrival, i.e., U^' n n . There is a metric tensor that allows to raise and lower index, 



(j)nn' _ (j) ( 3 3 \ _ f_iu-nm 

6 ~ e nn ~ \ n n' ) ~ [ L) 

With this metric we can now obtain the intertwiners for trivalent vertices with mixed ingoing/outgoing lines, 



(6) 
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n m i J2"i2 J3TO3 



(7) 
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-[ijimx j 2 m 2 , _ Jj 3 ) -\nini! j 2 m 2 j 3 m' 3 

v 33™3 ~~ fc m' 3 m 3 V 



(8) 



As an example of our conventions, let us consider, the evaluation of the Wilson net for the so called "theta diagranr 
spin network, 



+ 




We have therefore provided a definition of the Wilson net for any spin network. One simply constructs an invariant 
with holonomies and invariant tensors as the ones described above. 

With the above definition, the Wilson net has certain properties. In particular there is a dependence on the Wilson 
net on the orientation of the vertices. What we will now do is modify the definition in order to have an invariant that 
does not depend on the orientation of the vertices. This will correspond to the normalization of the Wilson net given 
by Witten and Martin. What we do is multiply the definition introduced up to now times a factor given by, 



V± = exp \± — \j x + j 2 + j 3 ]j </(2j 1 + l)(2j 2 + l)(2j 3 + l) (10) 

for each vertex in the spin network. With this factor, one can show that the Wilson net defined is invariant under 
changes from + to — type vertices. It is worthwhile noticing that in the work of Kauffman and Lins [ fl8| on invariants 
of spin networks the normalization is different. 



C. Knot invariants for spin networks 



We wish to study invariants that represent the transform of the Chern-Simons state into the spin network basis. 
That is, we are interested in expressions of the form, 



E(T,k)= J DAcxp(ikS CS )W(T,A) 



(11) 



where T is a spin network and W(T, ^4) is the Wilson net we introduced in the previous subsection. This kind of 
integral has been analyzed using monodromy techniques |lj,[l(| and variational techniques |19| , po|] . The result is a 
regular isotopic invariant of spin networks. The techniques do not give a unique answer for the invariant, but there 
are several possibilities, tantamount to having several definitions for the measure DA. Each possibility is uniquely 
characterized by prescribing a value for the invariant on the so called theta-net. The choice of Witten and Martin 
PLPI is, 



9(jl,j2,j3) 




\ 



\J AiA 2 A 3 



(12) 



with, 



Aj = E 




qi+ 2 — q o 



(13) 



and as we mentioned above, q = exp(27ri / k) . 
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The definition of the invariant E(T,k) is completely given in pOj ], one needs to supplement the above definitions 
with relations for recoupling and for crossings of the invariant. The resulting invariant is regular isotopic, meaning 
that it is not invariant under the elimination of twists. Specifically, 




= (_i)ii+A+Ja ex p (wr(/u +h 2 - ha)) E 



= exp (— 1-Kihj) E 



= exp (2Ttihj) E 




(14) 



(15) 



where hi = ■ 

These moves correspond to diffeomorphisms, therefore the invariant constructed is not invariant under diffeomor- 
phisms. In other words, the invariant is diffcomorphism invariant of ribbons and not of ordinary loops. One needs a 
framing (a prescription for assigning a ribbon to each loop) to have a well defined invariant. 

However, it was noted early on in the context of loops, and later in the context of spin networks ||, that the 
dependence on framing can be concentrated on an overall factor. In order to isolate this factor, one simply considers 
a power series expansion in k = (2m) /k and extracts the linear coefficient in k. That coefficient, exponentiated, is 
the overall factor. This construction is discussed in detail in ||], where it is shown that it indeed leads to an ambient 
isotopic invariant, in other words, a genuine diffeomorphism invariant function of loops. The resulting invariant ^] is 
then given by, 



E(T, k) = E(T, 0) exp (i> x (T)k) P(T, k) 



(16) 



where P(T, k) is a spin network generalization of the Jones polynomial (it reduces to it when T is a simple loop in 
the fundamental representation) and v\ is the first coefficient in the expansion in powers of k (for a single loop it 
reduces to the self-linking number). The factor E(T, 0), which corresponds to the evaluation of the invariant for k = 0, 
contains information about the coloring of the graph, and no information about the embedding. It can be thought of 
as the evaluation of the Wilson net for a flat connection. 

It was shown by Alvarez and Labastida |^1[ in the case of V being simple loops, and later generalized for links, that 
the invariant E(T, k) can be written as the exponential of a linear combination of primitive Vassiliev invariants, 



E(T, k) = E(T, 0) op ay {T)r l3 (G)k' 



(17) 



where ay(T) are the primitive Vassiliev invariants, which are only dependent on the embedding of the loop (they are 
independent of the gauge group of the Chern-Simons theory), di is the number of primitive Vassiliev invariants of 
order i, G is the gauge group (for our case SU(2)), and are group-dependent coefficients. We will assume that a 
similar structure appears for the case of generic SU(2) spin networks, specifically 



E(T, k) = E(T, 0) exp [ ^ v^T)^ 



(18) 



That is, we assume that the invariant is still given by the exponentiation of a set of "Vassiliev" invariants for spin 
nets, but we are not decomposing the expression in terms of "primitive" invariants, which is a structure that at present 
is not known in the spin network context. The Vi(T)'s are ambient isotopic for i > 1. 



1 Notice that since the definition of the theta-net involves square roots of the variable, the invariant is no longer a polynomial, 
as it was in the case of loops. It becomes a polynomial for the case of colored links, i.e. spin networks without intersections. 
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FIG. 1. The loop that defines the loop derivative 



By constructing these invariants of spin networks, we have bypassed one of the main obstacles that working with 
loops and links had in this context: how to comply with the Mandelstam identities. The invariants of spin nets we have 
constructed automatically take care of this. An important missing element in the spin network context is the idea of 
how "generic" is the set of Vi(T)'s. In the context of loops it is conjectured that the primitive Vassiliev invariants are 
enough to distinguish all knots. In the case of spin networks we are not working with primitive invariants and therefore 
it is questionable how generic the basis of invariants one is considering is. This is important if one is making the case 
that these invariants are the "arena" where one is going to discuss quantum gravity. If a decomposition in terms of 
primitive invariants of the exponential were known, the techniques we will develop will still be applicable. However, 
since up to present it is not known how to do this decomposition we will work with the v^s. It is worthwhile pointing 
out that the techniques we will introduce later in this paper to operate with loop derivatives and diffeomorphism 
constraints are geometrical in nature and not group dependent. Since it is known that all Vassiliev invariants can 
be constructed from the Chern-Simons integral with arbitrary groups, and our technique is not group dependent, we 
therefore have a de-facto method to operate on all Vassiliev invariants. For concreteness we will in the current paper 
concentrate on the case of SU(2). 



The loop derivative |j] is a differential operator in loop space that arises by considering two loops that differ by an 
infinitesimal element of area as "close" . It acts on basepointed objects (either loops or spin networks) by adding a 
path starting at the basepoint up to a point x, where it introduces an infinitesimal planar loop, and retraces back to 
the basepoint. The definition is, 



element of area spanned by the two vectors u and v defining the parallelogram one adds at the end of the path 7r. 
For a spin network, acting at a line on the network, the definition is exactly the same as for loops, except that the 
path becomes a line of the same valence as the line on which the basepoint lies, and the infinitesimal loop is also of 
the same valence as the path. 

We now wish to apply the loop derivative to the invariants we introduced in the previous section. A priori one 
expects that such a quantity does not exist. In particular, for a generic knot invariant, the loop derivative indeed 
does not exist. This is due to the fact that knot invariants are discontinuous functions in the space of loops. Loop 
derivatives change the topology of loops (for instance, they can remove intersections |ll]]), and therefore the limit 
defining the derivative is ill defined. What we will show here is that due to the properties of the invariants of Chern- 
Simons under deformations of the loops given by the skein relations, one can introduce a reasonable definition of the 
loop derivative for these kinds of invariants. It is similar to try to define the derivative of a discontinuous functions 
by allowing the derivative to take value in the distributions. We will analyze the consistency of this result with the 
properties of the invariants. The strategy is as follows. The invariants are defined as a functional integral of a Wilson 
net with a weight function. The only dependence on the spin net is in the Wilson net, so we will assume that the 
loop derivative of the invariant is equal to the functional integral of the action of the loop derivative on the Wilson 
net, with the appropriate weight factor, 



III. LOOP DIFFERENTIABILITY OF THE VASSILIEV INVARIANTS 
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FIG. 2. The loop derivative acting on a spin network 



Here we have assumed that the limit defining the loop derivative and the one defining the path integral are inter- 
changeable. 

The starting point of the calculation is the action of the loop derivative on a holonomy (in the j representation of 
SU(2)) associated to an edge of the spin network, containing the basepoint o, 



A o6 (<)^«(e^) m n = ^W(e^)^«(<)F«(x)y«(7r- i :)^)(ef)]- n , 
and the fundamental relation satisfied by the Chern-Simons state, 



47TZ S 
7 ^abc 7~\ 

k SA { c 3) (x) 



exp(ikS C s[A]). 



(21) 



(22) 



So the idea works exactly in the same way it did for ordinary loops |Tl],|2^] , the loop derivative acting on a holonomy 
produces an F a b, which can be re-expressed as a functional derivative acting on the exponential of the Chern-Simons 
form. The last step is to perform a formal integral by parts of the functional derivative and have it act on the 
holonomies of the spin network that were left by the action of the loop derivative. The end result is, 



A ab (n x )E(T, k) = -2k V e abc f dy c S 3 (x - y) x 



(23) 



< 



>cs 



where t/a are the SU (2) generators ^| in the j representation, and the expectation value is assumed to be taken with 
respect to the measure DAexp(iScs) and the dots refer to the fact that we just highlight the portion of the spin 
network where the loop derivative acts. It is understood that the products of holonomies continue until the net is 
closed and the appropriate traces are taken. A pictorial representation of the quantity within the expectation value 
is given in figure |2j. 

The above expression can be rearranged using the Fierz identity, 



(fc)Jsl'(j)Jr 



1=1 



j+k 
l=\j-k 



(2l + l) Pl (j,k) 



ks 



(24) 



jr 



where pi — 1(1 + 1) — j(j + 1) — k(k + 1), and where the dotted line denotes that the two vertices are at the same 
point, i.e. the line between the two vertices carries a holonomy equal to the identity. 
Using the rearrangement, we get, 



A a6 «)£(I» 



5^ K (2l + l) Pl (j,k)e abc fdy c S 3 (x-y)/ 

eger l=\j-k\ Je c \ J 




cs 



(25) 



2 Our convention for the SU(2) generators is to take the Pauli matrices divided by 2. This differs from other authors 
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and then using the basic recoupling property of SU(2) spin networks, 



P s 



we get, 



q r 



y = E(- 1 ) 2j ( 2m+1 ) JU 



ro=0 



F 



JP 



(26) 



j+k 2j 



A ob «)£(I» = - 



E K E E(-l) 2j '(2m+l)(2Z + lke ai)c f D dy c 6 3 (x-y)l 

eggr l=\j-k\m=0 Je c \ J 



B 




(27) 



cs 



and applying recoupling in the dotted circle in ( J27j ) and using the 6j symbols, one gets the final expression, 



A ab (^ )E(r, «) = E E I D d y cs3 ( x - y) x ^ k ) ( 




(28) 



cs 



where, 



A±(j,fe) = (-1) 



(_ ]_)'-(™=FW»)/2 

2=|j'-fc| 



(2m+I)(2Z + l)p; 



J J 
k k 



(29) 



So the end result is that the loop derivative acting on the invariant E(T, k) is nonvanishing only if the cndpoint of 
the path ir^ of the loop derivative falls upon one of the lines of the spin network. The end result is, up to a factor, 
proportional to the invariant E(V , k), where V is a new graph obtained by adding the path ir^ to the original graph 
r. Notice that the action of the loop derivative is covariant under diffeomorphisms, in the sense that a diffeomorphism 
would shift both the graph T and the path ir% and therefore would act as a diffeomorphism on the graph V . It is also 
worthwhile noticing that the general form of equation ( p8[ ) is true for any gauge group, the only differences would 
appear in the recoupling coefficients A and the corresponding irreducible representations associated with the graph. 
This emphasises the geometrical nature of the action of the loop derivative and may open possibilities for generalizing 
this construction for invariants that do not necessarily arise from SU (2) groups. The importance of this is that it 
appears that Vassiliev invariants for spin nets might arise as linear combinations of products of "primitive" invariants 
associated only with the topological embedding of the graph of the spin net, times some group-dependent factors 
that contain information about the valences of each line in the spin net. The loop derivative acts on such objects 
by ignoring the group prefactors and acting on the factors depending on the embedding of the spin net diagrams. 
In particular, if one acts on E(T, 0), since it does not have information about the embedding, the loop derivative 
automatically gives zero. 

The loop derivative as defined in this paper has several appealing properties that other differential operators in loop 
space do not have (see || for more details). One property we wish to emphasize is that the loop derivative satisfies 
Leibnitz' rule. It acts on a product of functions exactly like an ordinary derivative. This allows, when evaluating 
operators on products of knots (as for instance when we extract the frame dependent prefactor and the E(T, 0)), to 
perform explicit calculations. 

For instance, we can compute the action of the loop derivative on any Vassiliev invariant of type v n , by computing 
the logarithmic derivative of E(T, k). The final result is, 



2.; 



ex P I E M r m) - v r {T)) K r 



K=0 

(30) 



and the action is nonvanishing only if o, x fall on two different lines of T and the spin net T rn is obtained by taking T 
and adding to it the line tt^ with valence m. 
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IV. THE DIFFEOMORPHISM CONSTRAINT 



A. The operator 



In this section we will introduce the diffeomorphism constraint as a differential operator in loop space written in 
terms of the loop derivative. We will also show that the invariants we introduced in the previous section will be either 
annihilated by the constraint in the case of ambient isotopy invariants, or will transform appropriately in the case of 
regular isotopic invariants. 

Let us consider the diffeomorphism constraint of quantum gravity written in terms of Ashtekar's new variables, 

C{N)$[A] = J d i xN{x) a E b l {x)F l ab {x)^[A]. (31) 

Because this expression involves the product of two operators at the same point, in general wc need to regularize it, 
which we do via a point splitting function lim e ^o fe(x, y) = 8 3 (x, y), 

C(NMA] = limj d 3 x J d 3 2/ / e (x,y)iV Q ( 2 ;) J l-(x)[(7(^)^a fc (y)C/~ 1 K)]^[A] (32) 

where in order to preserve gauge invariance we have connected the F ab and E operators with holonomies along a path 
7r going from x to y. 

The above operator, when acting on a Wilson net, can be rewritten in terms of the loop derivative, as has been 
discussed in the context of loops in the past The explicit expression is, 

C(N)W(T,A) = lim f d 3 xMx,y) £ f d y b N a {y)A ba {^ y )W{Y,A). (33) 

One can explicitly check that the action of this operator is a diffeomorphism, provided that the connection A is 
smooth. In the limit in which the regulator is removed, the path 7r shrinks to a point and one ends with a loop with 
an infinitesimal closed loop attached at the point x. The addition of this closed loop is tantamount to displacing 
infinitcsimally the line of the original loop at the point x. 

We will now assume that the generator of diffeomorphisms has in general the form given by equation d33| ) in terms 
of the loop derivative and we will show that when acting with it on the invariants from Chern-Simons we get the 
correct result. This result is nontrivial, since in the path integral that defines the invariant there are contributions 
from distributional connections. 



B. Action on regular isotopic invariants 

To discuss in a cleaner fashion the action of the diffeomorphism operator on the invariants, we will consider diffeo- 
morphisms in which the vector N has compact support. This will allow us to focus on the action of diffeomorphisms on 
individual edges, on vertices, etc, each at a single time. It is immediate that a generic situation can be analyzed com- 
bining all results we will derive. Let us start with the action at an individual line, (we define C(N) — lim c _»o C e (N)) 



C e (N) 



= KitiJ) J d V b J dz c f e (z,y)e abc N a (z) 



(34) 



cs m=0 



z- 



' y + 



+ e(y-z) 



cs 



> 



cs 



where O are Heaviside functions that order points along the line of the spin net we are considering. One can easily 
define them introducing a parameterization. We now use a recoupling identity, 
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)■ 



cs 



(-l) 2i r . . , 

(2jTT) {j ' J ' m} 



(35) 



cs 



where {j,j, m} is the 3j symbol that is zero unless the j and m values satisfy the triangular relation at the vertices. 
We now remove the regulator, assuming the following regularization function, 



Evaluating the integral, one gets, 



fe(y,z) 



47T6 3 



Q(e-\z-y\). 



(36) 



2.1 



CS 



£: E x mU,f, 



Sit 



m—0 



(-i) 2j 

(2j + 1) 7 



6 



a6c " 



^(s) a 2/(s) c j/(s) 6 



|2)(^)I 3 



(37) 



cs 



where we have introduced a parameterization such that dy a = y(s) a ds and dots refer to total derivatives with respect 
to s. We can summarize the above result by saying that, for diffeomorphisms of compact support acting on lines of 
the spin net, we have that, 



C(N)E(T, k) = Kfx(J)w(N)E(T,K) 



(38) 



where fi(J) = Ylm=n j) (2j+T) ^ s a rec oupling factor, that depends on the weight of the line in question, and 

w{N) is the writhe introduced in the line by the action of the diffeomorphism along the vector N . We therefore see 
that the diffeomorphism has a well defined limit. The result is dependent on a background metric through the writhe, 
as expected for a regular istopic invariant. Notice that the result decomposes as a product of a factor depending 
entirely on the "coloring" of the spin net and another factor depending only on the embedding of the net in three 
dimensions. What we have therefore recovered is a formula that contains the information about the non-invariance of 
E(T, k) under the addition of twists, given by the skein relation (|l5|). This skein relation could be reobtained exactly 
by exponentiating the action of the loop derivative, as discussed in p2j, but we will not repeat the calculation here 
for brevity. 

Let us now consider the action of the diffeomorphism at an intersection. Here we get, in addition to the contributions 
we discussed for a single line, a contribution when the end of the path of the loop derivative intersects one of the lines 
entering the intersection. We have, 



C(N) 



3 2 Ja 

K E E KfJ^jp) / dy b / dz c f t (y,z)N a (z)e a 




cs 



(39) 



and using the recoupling identity, 




cs 



(_l)m+j-y-j a -j I 3 a Jf3 m 1 j 

j 30 JCL jj ' ' 



cs 



(40) 



and removing the regulator as before, we get, 



C e (N) 




2 3 

= K E E °m(ia,i/3,i 7 )ea6ciV a (0)^(0)^(0) 
CS a=l/3=ct + l 




CS 



(41) 
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where, 



Q m U,k,l)=Q(-ir +l - j - k \ I \ 7 \Kn(j,k) (42) 



with Q a constant given by, 



Q = ± r (sin 0- cos 0) (43) 

y a (0) | 2 cos 2 <j) + | y p (0) | 2 sin 2 - 2 1 Q (0) ■ y g (0) | cos sin 



4tt 7 



In the above equations we have assumed a parameterization with parameter equal to zero at the intersection. The 
geometrical meaning of this expression is that it vanishes if the volume spanned by the two tangents at the intersection 
and the derivative of the vector N along one of the lines entering the intersection vanishes. If this volume is nonzero, 
we start to reconstruct infinitesimally the twist at the intersection covered by the skein relation (Q). Again, one 
could recover it exactly exponentiating the action of the loop derivative as described in p2[ . 

C. Action on ambient isotopic invariants 

Here we will show that the diffeomorphism operator we introduced vanishes when acting on the ambient isotopic 
invariants, concretely P(T, k). Therefore, it annihilates every Vassiliev invariant. In order to compute the action, let 
us recall equation (|l8|), 

E(T, k) = E(T, 0)(1 + vi(T)n + {2v 2 (T) + t>i(r) 2 )y + • • •) (44) 
and analyze the action of the diffeomorphism on E(T. k) order by order in k. To order zero we get, 

C(N)E(T, 0) = (45) 

which is correct, since E(T, 0) does not depend on the embedding of the spin net, just on its coloring, and it is 
immediately annihilated by the loop derivative. At the next order, we have, 

C(N) Vl (T) = fi{j)w(N) (46) 

where we have assumed the diffeomorphism as acting on an isolated line of valence j, a similar formula holds for the 
intersections with a different coloring weight. As before w(N) is the writhe induced in the line by the vector field N. 
If we now take advantage of the fact that the diffeomorphism operator satisfies Leibnitz' rule, we see that combining 
(p8|) and (Em we get, 



C(N) expi-KV! (j))E(T, k) = 0, (47) 
and therefore we see that indeed the diffeomorphism constraint annihilates all Vassiliev invariants, namely, 

C(N)P(T, k) = C(N) exp (j^ = «• ( 4 «) 

Since the diffeomorphism constraint annihilates all Vassiliev invariants, it also annihilates any function of the 
Vassiliev invariants. As we mentioned before, there are indications that Vassiliev invariants may constitute a basis of 
diffeomorphism invariant functions of loops. What we have shown here is that the definition of the loop derivative we 
introduced in the spin network context is compatible with that fact, namely it naturally leads to a diffeomorphism 
constraint that annihilates explicitly all Vassiliev invariants. 
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V. TOWARDS A HAMILTONIAN CONSTRAINT 



Let us now consider the Hamiltonian constraint of quantum gravity, possibly with a cosmological constant A, 



H(M) 



d A xM 



e tjk E, E,F k h 



e^E^E^ 



(49) 



This expression corresponds to the "doubly densitized" operator. We chose to do this since we have experience with 
this operator from the context of loops || and it is easily expressed in terms of the loop derivative. We will later 
suggest how one could transform this expression into a single densitized version. We need to regularize the operator, 
which again we do through a point splitting. There are many options as to how to exactly point split. The way we 
do it, for the A = piece, is, 



. («, v, w, x)H(u, v, w) = 9e , e .,Au, v, w, x)e^ k E t (u)EJv)[U(^)F ab (w)U(Trl)] k . 



(50) 



So we have chosen to join the F ab with one of the £7's via a holonomy. There are many other possibilities (for instance, 
joining all operators via holonomies, which would yield a gauge invariant regularization) and they all yield the same 
classical expression in the the limit e — » 0. A possible regulating function is defined as, 



f(u,V,W,x) = f € (x,u)f e (x, 



—)fe(X, — 



- w 



(51) 



with / defined as in the diffcomorphism constraint. This regulating function is the same as the one usually consid- 
ered for the "Ashtekar-Lewandowski" volume operator, fl],|23||. It is worthwhile pointing out that for the following 
discussion, this issue is not important, since we will discuss results that do not require to remove the regulator. 

To determine the action of the Hamiltonian in terms of spin networks, we consider the action of the operator we have 
just defined on a Wilson net. As is well known, this operator only acts at intersections of the net. At regular points 
it gives rise to "acceleration" terms, which we will omit since they vanish on diffeomorphism invariant states [p4],p[. 
The nonvanishing action comes from the triads acting at two points on different strands i, k entering an intersection 
V, 



H(u, v, w)W 




i+k 



= ~2 2 E (2m + l)p m (i,k) dy a S 3 (u~y) dz b 6 3 (v-z)x (52) 

&keV m=\i-k\ ei eh 




\ 



J \m 



/ 



where in the last diagram all the structure occurs "at a point" and points yi are identified. 

Having the Hamiltonian in terms of the loop derivative allows to evaluate its action on any of the invariants we 
have been discussing. Let us start by showing that the Hamiltonian constraint with a cosmological constant term 
annihilates the invariant E(T, k). In order to do this, we need to derive an expression for the operator detq, q being 
the spatial metric in terms of spin nets. The construction completely parallels that of the Hamiltonian constraint, 
the operator detq is obtained by replacing F* b in the Hamiltonian with e a bcEf. In order to be possible to have a 
cancellation with the Hamiltonian, one has to regularize both operators in a similar manner. So we choose, 



flr e>e / >e » (it, v, w, 2c)det(u, v, w) = ^g^\e"(u, v, w, x)e lJk E i (u)Ej(v)[U «)^ c (w)[/«)] fc e Qbc 



(53) 



The point-splitting (i.e. where one decides to insert each holonomy) could have been made in a different way (e.g. 
PIP]), but one can show — using recoupling theory — that inserting additional holonomies leads in the limit when one 
removes the regulator to the same result for a given regulating function. 
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The end result for the operator is, 

/ , \ 

detg(u, v, w) 1 ^ 




be X 



(54) 



i,k,l m,n 



\ J i h) 



dy a 5 i {y,u) I dz b S%z,v) I dt c 8 6 {t,w) x 

( \ / \ ( 

- * 





The vertex on which the determinant acts can be of any order, if it is higher than trivalent, the determinant is a 
sum of the above contribution evaluated for all possible triplets of lines. All linkages in the figure are assumed to 
be happening "at a point" , so recoupling theory can be applied. An interesting observation is that if the vertex is 
trivalent, simple recoupling identities show that both contributions are equal and therefore the determinant of the 
metric vanishes on such vertices (as one would have suspected, given that the volume operator vanishes too). 
Using the definition of the above operators, one can start obtaining results. The first thing we notice is that, 



ge,e',e" (u, v, w, x)H(u, v, w)E(T, k) = YIkAqh^i i£ » (u, v, w, x)detg(u, v, w)E(T, k) 



(55) 



and we therefore see that the generalization of the Kauffman bracket to spin networks E(T, n) is a solution of the 
Hamiltonian constraint of quantum gravity with a cosmological constant. This is straightforward to see by considering 
the action of (|52] ) and making explicit the action of the loop derivative, which adds a line to the diagram in such a 
way that one ends with the same diagrams as in (|54|) . 

Another interesting result is to notice that the Hamiltonian vanishes when acting on Vassiliev invariants at trivalent 
vertices. This can be straightforwardly seen by recalling the action of the loop derivative on a Vassiliev (|3(]), and 
noting that at a triple intersection, using recoupling, v p (Tj) = v p (T). We are allowed to use recoupling because 
the loop derivative acts by adding a line "at a point" when evaluated in the action of the Hamiltonian constraint. 
Therefore the Hamiltonian constraint vanishes at trivalent vertices. 

Another important result concerning the action of the Hamiltonian on Vassiliev invariants can again be concluded 
from studying the explicit action of the loop derivative (|30|). Combining this expression with (52) one immediately 
concludes that the action of the "doubly-densitized" Hamiltonian we are considering here is proportional to the 
"regularized volume" (strictly speaking it is a volume squared) spanned by three of the tangents entering at an 
intersection |2al, 



ijk 



Vol reK (a;,e, e ,e ) = e abc / dy a / dz a / dw a g e ^^>(x,y, z,w) 



Concretely, acting on a four-valent vertex labelled by J, the action of the Hamiltonian is, 



(56) 



> (x, u, v, w)H douh i c (u, v, w)v n ( 




E E vo, 

/ i,j,fe£Vertex 



ijk 



(x,e,e',e")C IJ V n 




(57) 



with C IJ is a proportionality factor depending on the valences of the strands entering the intersection. Similarly, 
the ranges of the sum in I is determined by the type of intersection. The invariant V n -\ is a non-primitive Vassiliev 
invariant of order n — 1, it will involve sums of products of primitive Vassiliev invariants of lower orders. 

This proportionality opens the attractive possibility of defining a "singly densitized" Hamiltonian by "dividing 
by the regularized volume" . Having a singly densitized Hamiltonian can potentially lead to a much better defined 
operator, since its action could possibly be cast in a background-independent manner, as proportional to a Dirac 
delta. If this construction works, it might therefore be possible to cast the action of the Hamiltonian constraint in 
terms of an expression along the lines of, 
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H, 



single 




S(x - Vertex) ^ d /J K-i 




(58) 



with d IJ another constant factor. Having an operator with such a simple action is not only remarkable, but will open 
the possibility of further consistency checks, like computing the constraint algebra. All these issues are currently under 
study and evidently further work is needed to complete this program. The computation can be done "on shell" by 
showing that the Hamiltonian commutes on Vassiliev invariants using the expressions introduced above. It could also 
be done "off shell" since one has an expression for the Hamiltonian acting on any function that is loop differentiable, 
and therefore verify commutation relations with the diffeomorphism constraint. 



VI. CONCLUSIONS 



The doubly-densitizcd Rovelli-Smolin Hamiltonian has been studied, written in terms of the loop derivative, since 
the early days of the new variable formulation in terms of loops. The early work suffered from three main drawbacks: 
a) the loop derivative was not a well defined operator on diffeomorphism invariant states; b) the quantum states 
were difficult to construct and the action of the constraint to check explicitly due to the presence of the Mandelstam 
identities and c) the Hamiltonian was regularization dependent, potentially having problems at the time of computing 
the quantum commutator algebra of constraints. In this work, largely by going to the language of spin networks, 
we have made progress in all three points. First of all, we were able to introduce a loop derivative that is well 
defined when acting on Vassiliev invariants. The loop derivative leads, while input into the usual definition, to a 
well defined diffeomorphism constraint, with exactly the kind of action on Vassiliev invariants and its regular- isotopic 
counterparts, as one would expect. The presence of this reasonable-looking loop derivative operator allows us for 
the first time to actually have a version of the Rovelli-Smolin Hamiltonian written in terms of the loop derivative 
that actually is a well defined operator on knots. We have checked that several of the formal results obtained in the 
loop language still go through with the new operator. Better yet, since all the results are in terms of spin network 
states, we need not concern ourselves anymore with the Mandelstam identities. Finally, since the Hamiltonian has 
a relatively simple action on Vassiliev invariants, closely related to the volume operator, the possibility is opened of 
defining a single-densitized Hamiltonian, just by dividing by the operator representing the determinant of the metric. 
This latter possibility still needs to be studied in detail. However, preliminary results appear as promising. The single 
densitized operator appears to have a very natural action on the space of knots. This opens the hope that a consistent 
set of constraints, satisfying the appropriate commutator algebra, might be found. One could think of computing the 
commutator algebra off-shell, very much along the same lines as in the formal computations of ||, except that now 
the operators involved are also well defined on-shell (i.e. when the wavefunctions arc knot invariants) . The fact that 
all this occurs while working on the space of Vassiliev knot invariants based on spin networks suggests that this might 
be an appropriate setting for the discussion of canonical quantum gravity. 
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